Non-Markovian Quantum Trajectories of Many-Body Quantum Open Systems 



Jun Jingi*, Xinyu Zhao\ J. Q. You^t, Ting Yu^* 
^ Center for Controlled Quantum Systems and Department of Physics and Engineering Physics, 
Stevens Institute of Technology, Hoboken, New Jersey 07030, USA 
^Department of Physics and State Key Laboratory of Surface Physics, Fudan University, Shanghai 200433, China 

(Dated: December 3, 2010) 

We have derived for the first time a set of non-Markovian quantum trajectory equations for many- 
body and multi-state quantum open systems. The exact time-local (time-convolutionless) quantum 
trajectory equations can be used to simulate quantum dynamics of coupled open systems interacting 
with a bosonic environment at zero temperature. Our general results are explained and illustrated 
with several atomic and optical open systems. 
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Introduction. — Dynamical aspect of quantum open 
systems has been investigated for many years and devel- 
oped in different formalisms |TH9]. Typically, the state 
of an open system is described by a density operator 
governed by a master equation that plays a pivotal role 
in the fields of quantum optics and quantum informa- 
tion. When the environment is a structured medium or 
the system-environment interaction is strong, Lindblad 
Markov master equations are prone to fail badly. Then, 
non-Markovian master equations or their alternative non- 
Markovian quantum trajectory approaches must be used 
[TUHT3]. Despite extensive efforts, deriving a master 
equation without invoking the Born-Markov approxima- 
tions (weak couplings and mcmorylcss environment) has 
achieved only limited success in practice, and is mostly 
restricted in a few isolated models such as the quantum 
Brownian motion model J^, a leaky cavity model [15, 16J, 
a phase damping model [T7j and a two-level atom coupled 
to a quantized radiation field '18] etc. For a many-body 
open system such as an iV-qubit system or a multi-state 
system coupled to a fully quantized environment, the ex- 
istence of non-Markovian dynamic equations [191 such as 
master equations is still largely unknown. 

It has been shown that the quantum state diffusion 
(QSD) equations for the pure states can formally solve a 
quantum open system coupled to a bosonic environment 
irrespective of environmental memory, coupling strength 
and the spectral density [5DH21j. The generality of the 
QSD equation has given it appeal as both a numerical 
and theoretical tool for a non-Markovian open system. 
On the other hand, the time-nonlocal feature of the non- 
Markovian QSD equations has been a major obstacle in 
general implementations of the trajectory formalism for a 
realistic physical system. Therefore, it is highly desirable 
to develop a time-local non-Markovian formalism that is 
applicable to many-body open systems. 
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In this Letter, we report explicit constructions of ex- 
act time-local QSD equations for several many-body and 
multi-state quantum systems. Moreover, we show that 
these exact time-local QSD equations contain only poly- 
nomial noise terms. To the best of our knowledge, these 
time-local QSD equations for many-body open systems 
have never been derived before. Such non-Markovian 
quantum trajectory equations are capable of describing 
the transition from non-Markovian to Markov regimes for 
a genuine many-body quantum system. 

The model and the quantum state diffusion equation. — 
A generic quantum open system in the system plus envi- 
ronment framework can be be written as (setting h — 1): 

Htot — Hgys + Hint + Hcnv (1) 

where -ffgys is the Hamiltonian of the system of interest. 
Hint = X]k(5k^cik dk^^'^k) is the interaction Hamil- 
tonian and -ffcnv = X^k'^kojj.ak describes a quantized 
field (environment). Note that L is a system operator 
characterizing the mutual interaction between the sys- 
tem and the environment. At zero temperature T = 0, 
the environmental correlation function is determined by 
the noise operator B{t) = X^k ^kOke"'^''* in the interact- 
ing picture: a{t,s) = {0\[B{t) + BHt)][B{s)+B^s)]\0) = 
Ekl5kpe— (*-^). 

Here we use |4'tot(0) to represent the state of the total 
system at time t, then the reduced density operator pt for 
the system of interest is given by pt = Trenv[|^tot>(*tot|] 
obtained by tracing over the environmental degrees of 
freedom. If the system and its environment are initially 
uncorrelated, it has been shown that the density opera- 
tor for the open system can be decomposed as a set of 
continuous quantum trajectories living in the system's 
Hilbert space, denoted by ipt{z). The trajectory ipt{z) 
is governed by a linear stochastic Schrodinger equation, 
termed linear QSD equation [20l [21] : 

d S 

-Mz) = {~iHsys+Lzt-L^ / dsa{t,s) — )Mz) (2) 
at Jq ozs 

where zt = —i X^k Sk^^^e*"''* is a complex Gaussian pro- 
cess satisfying M[zt] = M[ztZs] = 0, and M[ZfZs] = 
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a{t, s). Here M[-] denotes the ensemble average over the 
classical noise Zt- Note that a{t,s) is the correlation 
function defined above. The quantum trajectory tptit) 
is designed to recover the density operator of the system 
by taking ensemble average: pt — M[\'tpt{z)){ipt{z)\] — 

The Diosi-Gisin-Strunz QSD equation ([2]) is a remark- 
able result since it is completely general irrespective of 
spectral density of environment and coupling strength. 
The appearance of functional derivative in Eq. ^ is 
strongly reminiscent of the convolution kernel appear- 
ing in the Nakajima-Zwanzig master equations obtained 
from projection operator technique [25]. Thus, for the 
purpose of practical applications, it remains unresolved 
whether one can recast the existing functional derivative 
into a time-local form. 

The ansatz to replace the functional derivative 
S'ipt{z)/6zs = 0{t, s, z)iptiz) combining with the consis- 
tency condition jf-^^ = yields the equation of 
motion of the 0-operator piT - lM] : 

^^Hi/sys + L.*-i^O,0]-Ltg. (3) 

where 0{t, z) = J* a{t, s)0{t, s, z)ds and the initial con- 
dition 0(s, s, z) = L is satisfied. Once the 0-operator 
can be explicitly constructed, then the original non-local 
QSD equation ([2| takes a time-local form: 

^Vt(^) - {-^H^ys + Lzt ~ L^O{t, z)) iPtiz)- 

Clearly, finding the solutions to the nonlinear operator 
equation ^ is by no means trivial. Up to now, only a few 
isolated examples have been explicitly solved [211 [521 [IS] ■ 
Here, we have derived the time-local QSD equations for 
a large class of many-body and multi-state systems that 
are of importance in quantum optics, atomic many-body 
physics and quantum information science [27j . 

For numerical simulations, one must use the nor- 
malized pure states ^/'((z) — where = 
{'ipt{z)\tl)t{z)) , which is governed by a nonlinear version 
of the QSD equation: 

- At{L^)0{t, 5)^Pt + {At{L^)0{t, ~z))ti^u 

where Ai(A) = A — {A)t for any operator A, {A)t = 
(V'tl^lV't) denotes the quantum average, and zt = zt + 
/p a*{t, s){L^)sds is the shift noise. For all the numerical 
results to be presented below, we alway use the normal- 
ized nonlinear QSD equation. 

Time-Local Non-Markovian QSD equations. — 
Equation ([S]) manifests that the 0-operator may be 
determined by the commutation relations involving Hgys 
and L together with a set of basis operators. 



In general, the 0-operator may be expanded as [25] . 

oo 

0(i,s,z) = 0(")(t,s) + ^OW(i,s,z) (4) 

k=l 

where 0^°'>{t,s) ^ J2j fji^^ s)of is the noise-free 
term; 0(i)(i, s,z) = J2j lo p'j^H^^ Si)zs^dsiO'j^^ is 
the linear-noise term, and in general O^'^^' (i, s, z) — 

^3 /q ■ ■ ■ Slpf\^^ s, si, • • • , Sfc)zsi ■ ■ ■ Zs^dsi-dskOf' 
contains fcth-order polynomial noises. Note that all the 
basis operators O^^^ and functions fj^P'!'^ are noise-free. 

(k) 

Moreover, the basis operators are time-independent. 

The equations for fjjpf^^ may be obtained from Eq. 
Below, we shall show that the polynomial 0-operators 
can be explicitly constructed for several many-body 
models ranging from discrete to continuous variable 
systems. 

Theorem 1. (A^-cavity model) For the open system 
model (1) with Hgys — + -|- 

(-^ + 1 — 1) and L — X^jLi then the exact 
0-operator is noise-independent. Explicitly, it is given by 
0{t,s,z) = 0'-'^'(i, s) = '^f^i fj{t, s)aj where fj satisfy 
a set of differential equations I28j . 

The proof and explicit construction of the 0-operator 
for iV-cavity model are straightforward (see [2H])- Here, 
we provide an example of the 0-operator with TV = 
3: 0{t,s,z) = fi{t,s)ai + f2{t,s)a2 + f3{t,s)a3, 
where fj{t,s), [j = 1,2,3) satisfy the following 
equations: §^fj{t,s) = iujjfj{t,s) + Aj_i/j-_i(i, s) -I- 
iXjfj+i{t,s) + fkiti s)Fj(t) with the notation 

Fj{t) = dsa{t, s)fj{t, s). The initial conditions are 
given by /i(s,s) = /2(s,s) = /3(s,s) = 1. 

Theorem 2. (Angular momentum model or iV-level 
model) For the open system model (1) with Hgys — (^Jz 
and L — where Jz and J_ are angular momentum 
operators satisfying [Jz,J±] ~ ±J±, [J_,J±] — +Jz, 
then the exact 0-operator 0(t, s, z) can be explicitly con- 
structed. More precisely, if I is the angular momentum 
quantum number (integer or half-integer), then the O- 
operator expressed in |^|) has the basis operators O^^'^ — 
J^- with k ^ 0, ■■■ ,21 - 1 and j = 1,2, ■■■ ,21 - k. 

Let us emphasize that it is possible to construct the 
explicit 0-operators for the multi-state systems with ar- 
bitrary energy spacings. The details of the proof for The- 
orem 2 can be found in [5H] . It is instructive to provide a 
simple example of the 0-operator of a spin-2 system (or 
a five-level system). Explicitly, the 0-operator may be 
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written as: 

0{t,s,z)= 

+ E f Pf\t,s,s,)z,,ds,0f^ 

+ Y pf\t,S,Si,S2)Zs^Zs2dsids20f^ 
= 1,2-' •'° 

+ jjj pf\t,s,si,S2,s:i)zs^Zs^Zs^dsids2ds3,of\ 

(k) 

where the time-independent basis operators Oj are 
given by of = J^^J-, (j = l,---,4); Of^ - 
Jr'JlAj = 1,2,3), Of^ = Ji-'Jl,U - 1,2) and 
= J1. The differential equations governing the co- 
efficients fj , p^'^'' are explicitly given in [28] . 
Theorem 3. (A^-qubit model) For the open system 

model (1) with H,ys = |Ef ^^^^ and L = ^ 
then the exact O- operator can be explicitly constructed, 
and contains (N — l)-order noises. That is, of ^ =0 if 

k;^ N. 

The general proof of Theorem 3 for an arbitrary N 
is rather cumbersome. Here we have proved Theorem 3 
for iV ^ 7 by explicit construction of the 0-operators. 
If we use m{N,k) (0 ^ k ^ iV — 1) to denote the 
number of terms with fc-fold noise integration in the 
0-operator, we have the following relations shown in 
Table |lj m{N,k) = m{N - l,fc - 1) and m(7V, 0) = 
m(7V - 2, 0) + iV with m(l,0) = 1. Table |l] lists the num- 
bers m{k,N) for seven cases. Obviously, when iV = 1, 
the 0-operator just reduces to 0{t, s, z) — f{t, s)cr_ 
for a single qubit case [HJ [53] with Hgys — '^Uz and 
L — a-. For N — 2, it is easy to check that the exact O- 
operator for the two coupled-qubit model has five terms: 
Oit,s,z) = E^ti/j(i,5)0f -t-i/oP(i,5,5i)z,,dsiOl''. 
where O^f^ = a^, ^ a^, 0^°) = a^a^, ^ 

a^af , and o[^'' = a^a^. As an application, we point 
out that the exact two-qubit time-local QSD equation 
allows us to calculate the entanglement evolution of the 
density matrix constructed from quantum trajectories 
[55]. Similarly, the exact three-qubit 0-operator can be 
constructed explicitly (See [5H]). 

Our result in Theorem 3 regarding 0-operators can 
readily be generalized to coupled A^-qubit systems, which 
is expected to generate quantum coherence of qubit sys- 
tems, including entanglement generation. The deriva- 
tion is also applicable to the iV-qubit systems where each 
qubit has a different frequency i?sys = J2j i-^jci'''' /2. Fur- 
thermore, iV-qutrit or iV-qudit models can be solved in 
a similar way. 

Numerical results and discussions. — Below, numeri- 
cal results for the non-Markovian quantum dynamics of 
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TABLE I: This table summarizes the explicit constructions 
of the 0-operator. The notation fe stands for the orders of 
noises contained in the 0-operator for an A'^-qubit dissipation 
model (iV = 1,2, ••• ,7). For = 2, m(2,0) = 2,m(2,l) = 
l,m(2,fe) = 0, (fe > 2). 

many-body systems are presented. For simplicity, and 
for recovering the Markov limit, we assumed that the cor- 
relation property of the noisy environment is described 
by the Ornstein-Uhlenbeck process: a{t,s) = ^e"'*''*"*'. 
Clearly, the Ornstein-Uhlenbeck noise approaches to the 
Markov limit when 7 — cx). It should be noted that our 
exact time-local QSD equations are valid for arbitrary 
types of correlation functions. 




FIG. 1; (Color online) Non-Markovian quantum interference 
of a four-level system (spin-3/2) with //ays = <^Jz and L — J-. 
Initially, the system is in the state ipo = | (| 1) -I- 12) -I- 13) -I- 14)), 
where |j)'s, j = 1,2,3,4 are the different levels for this four- 
level system. The time evolutions of coherences |p23| and 
P14I and population terms pn and P44 are calculated over 
1000 trajectories. The blue solid curve stands for 7 — 0.3, 
the red dashed curve for j — 1.0 and the black dot-dashed 
curve for 7 — 3.0. 

The first example provides a simple picture of the effect 
of the non-Markovian atom-cavity interaction on quan- 
tum interference of a four-level atomic system. Both co- 
herence decay and the population transition times from 
a higher level to a lower level are significantly different 
from the Markov limit as shown in Fig. [T] Our plots 
show that, in a non-Markovian regime, both coherence 
(Fig. la, lb) and population (Fig. lc,ld) exhibit decay 
and revival patterns before eventually decaying to their 
final values. So the loss of coherence has effectively been 
suppressed due to the memory effect of the environment. 
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(a) F=0 



(b) Y=0.3 




and atomic dynamics in photonic crystals. Furthermore, 
it would be of interest to extend the current methods to 
the case of hybrid quantum systems consisting of both 
continuous and discrete variables. 
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FIG. 2: (Color online) Evaluating fidelity under the influence 
of non-Markovian amplitude damping channel via quantum 
trajectories (1000 realizations). The initial state is the Werner 
state with parameter F: W — ^Ig + (1 — F)\ip){tp\, where 
IV) = l/\/3(1100) + |010) + i001)). (a) Fixed initial pure state 
with F = for different memory time r = I/7; (b) Fidelity 
for different mixed states for the same channel with 7 = 0.3. 

The second example shows how a three-qubit system 
evolves when coupled to a common multiple-mode en- 
vironment. The exact time- local QSD equation can be 
derived explicitly (For details, see jH]). Our numerical 
simulations with the zeroth-order 0^''-'(t, s) and the first- 
order terms 0^-^\t, s, z) reveal the novel features of co- 
herence dynamics measured by quantum fidelity. Fig. 3a 
shows the plot of hdelity against time for different en- 
vironmental memory times r = 1/7,7 = 0.3,1.0,3.0. 
Clearly, the shorter memory times give rise to better 
fidelity at time t. The result suggests a rather inter- 
esting feature that quantum coherence can typically sur- 
vive longer in a non-Markovian dissipation environment. 
Another interesting feature arising from this four-level 
system is that, for a strong non-Markovian environment 
with 7 = 0.3, the quantum fidelity is closely related to 
the degree of entanglement of initial three-qubit Werner 
states measured by the parameter F (Fig. 2b). As the ini- 
tial state approaches separable state, the state becomes 
less affected by the environment reflecting the fragile fea- 
ture of an entangled state [3D] . 

Conclusion and outlook. — We have established the ex- 
act non-Markovian QSD equations for several important 
many-body and multi-state quantum systems including 
A^-qubit systems and A^-level atomic systems coupled to 
a common bosonic environment. The results can be ex- 
tended to interacting qubits and high dimensional sys- 
tems. Our findings will have many applications in many- 
body atomic coherence dynamics and quantum informa- 
tion science as illustrated by two examples in this Letter. 

Our results can motivate other lines of research. 
Clearly, it is important to apply the exact many-body 
QSD equations to the important non-Markovian physical 
systems such as an atomic ensemble in an optical cavity 
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